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Abstract 

We extend the Majumdar-Papapetrou (MP) solution of the Einstein-Maxwell (EM) equations 
which is implied generally for static electric charge in non-rotating metrics to encompass equally 
well magnetic charges. In the absence of Higgs and non-Abelian gauge fields, 'dyonic' is to be 
understood in this simpler sense. Cosmologically this may have far-reaching consequences, to 
the extent that existence of multi-magnetic monopole black holes may become a reality in our 
universe. Infalling charged particle geodesies may reveal, through particular integrals, their inner 
secrets which are screened from our observation. 
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I. INTRODUCTION 



A relatively simple, yet interesting class of Einstein-Maxwell (EM) solution was given 
long ago by Majumdar jlj and Papapetrou [2], which attracted attention in various contexts 
including that of multi-black holes. The isotropic form of the line element with all inclusive 
metric function, which determines also the static electric potential makes this solution unique 
among EM solutions known to date. The metric function that generates the space time 
satisfies the Laplace equation. The linearity of this latter equation leads automatically to 
the multi-center solutions at equal ease. Each centre satisfies all the requirements necessary 
for black holes and as a matter of fact the multi-centre solution can be interpreted as a 
multi-black hole solution - which is otherwise extremely difficult to obtain analytically. For 
this accomplishment we are indebted to the Majumdar-Papapetrou (MP) form of the metric 
[3)]. A significant extension of the MP solution was to include time dependence through a 
cosmo logical constant J4]. This latter form of the metric paved the way toward black hole / 
brane collisions in higher dimensions {5, 6]. 

In this paper we wish to contribute to the MP solution by adding magnetic charge along- 
side with the electric charge. To our knowledge, MP space time has been considered so far 
only with a static electric field described by the potential = 5* A, in a diagonal metric. 
It is known that inclusion of rotation creates natural magnetic fields from the static electric 
charges [7|- Yet, by remaining in the static, non-rotating metric and adding a magnetic 
charge to the electric charge - i.e., a dyon - seems escaped from attentions. Let us note that 
this should not be confused with solutions such as Reissner-Nordstrom (RN) in which mag- 
netic and electric charges are treated on equal footing. Our magnetic charge lacks spherical 
symmetry since one of the axis (i.e. the z-axis) is singled out which is more apt for multiple 
axial superposition. The - dyonic-black holes consist of both electric and magnetic charges 
coupled together. Radial geodesic analysis of electrically charged particles infalling such 
black holes will exhibit different behaviors to aid in detection of such magnetic black holes. 
Under the light of such magnetic objects the natural question arises: Do such black holes 
serve as the storage of magnetic monopoles which are elusive in our observable universe? 

Although it is a matter of formality to extend our results to arbitrary number of black 
holes and to higher dimensions, we shall restrict ourselves to 4— dimensions and consider 
the example of 2— centre black hole as an example. The problem of horizon smoothness for 
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multi-black holes, and the issue of stability are two of the problems that we shall not address 
in this paper. 

Organization of the paper is as follows. In Sec. II we present the solution of EM equations 
with both electric and magnetic charges. Geodesic analysis follows for particular boundary 
/initial conditions. We extend our discussion, through perturbation, to the case of 2-black 
holes located along the z-axis. We complete the paper with Conclusion in Sec. IV. 



II. INTEGRATION OF THE EINSTEIN-MAXWELL EQUATIONS 



0,6 



We start with the Majumdar-Papapetrou line element in 4— dimensions [l|, |2| given by 

ds 2 = --^dt 2 + tt 2 (dx 2 + dy 2 + dz 2 ) (1) 
in which Q is a function of x, y, z and t. Our electromagnetic multi-centre potential ansatz 



is 



a = Ut + E r , Pi [ z 2 , f — 72T [(* - x i) dy-(y- yd dx] (2) 

where 

n = \j{x- Xi) 2 + (y- yij 2 + (z- z^ 2 , 

e is a constant such that < e < 1 and Pi stands for the magnetic charge of the i th black 
hole. The electric charge Qi of the i th black hole will be defined below (Eq. 22). The 
electromagnetic field two-form is given by 

F = dA = e i^dtdx + ^dtdy + ^dtdzj +E^" l( x ~ x i) dydz + (y - jji) dzdx + (z - z { ) dxdy] 

(3) 

with its dual 

*F = e (Q x dydz + Q y dzdx + Q z dxdy) + ^2 3 ^ 2 [(x — x^ dtdx + (y — yi) dtdy + (z — Zj) dtdz] 

(4) 

in which Q x , Q y , ... denote partial derivatives and dx % dx^ implies wedge product. Concerning 
Maxwell's equations, we have 

d{*F)=0 (5) 
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leading to 
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where Q xt = and so on. Eq. (6) is the usual Laplace equation whose simplest solution 
can be written as 

n = w(*) + £- (7) 

i T'i 

in which oj (t) is a function of time and Cj are constants to be identified. A substitution into 
the rest of Maxwell equations implies 



E3(y-w) £3(*-*) 

which is easily satisfied provided C* = \Pi for a constant A. From (3) and (7) we find the 
field tensor F^ u and the energy momentum-tensor as 

t; = 2F^\F uX - ±fs; (9) 



in which 



and 



F = F m F^ = -2 (Fg + F% + Fl) + A (F 2 xy + F 2 XZ + F 2 yz ) , (10) 



Px V 2^ 3 5 ^zz 2^ 3 5 ^J/z 2^ 3 \ lL ) 

i 1 i i 1 i i 1 i 

TP - ^ x TP - TP - ^ z 

The non-zero components of and G 1 ^ are tabulated in Appendix la and lb, respectively. 
A solution to the Maxwell Equations (5) follows once we set (8) to a constant (say A) i.e., 



or equivalently 



E Pdx 3 — = (13) 



r 



~^r~ = T' (15) 



Consequently, the field tensor components become 



F = 5i F = _^hi p = ^ (iR) 

1 xy ^ i 1 xz ^ 5 yz ^ j l xu / 



t-, _ jp _ ^z 



and as a result the energy-momentum components take the form given in Appendix lc. 

One may also substitute (from the Appendix) into the tt component of the Einstein's 
equation with the cosmological constant A to obtain 

Tl = G\ + A ->• (17) 
This is satisfied if we make the choices 

e 2 + ^ = 1 (18) 

and 

3Q 2 t = A. (19) 
The latter equation implies (from (7)) that 

u>(t) = ±^t + C (20) 

with an integration constant Cq that is disposable with the choice of origin of time. The 
rest of the Einstein's equations turn out to be satisfied all, by virtue of (18) and (19). In 
conclusion, we obtain the solution as 

n = ± A +Q)+? (21) 



which clearly for r = r, we have the location of the i th black hole with effective charge and 
mass equal to |PjA| . With reference to our potential ansatz (2), we observe that for e = we 
have the pure magnetic charge Pj = (mass) . To define the electric charge Qi we integrate 
the Maxwell equation in accordance with 

j E i .dA i = AitQ i (22) 

over the i th sphere to get Qi = ^^_ t2 Pi- It is clear from this definition that for e = 1 we 
must take p — > to have a meaningful electric charge, this is indeed the case as given in 
the sequel. For a single black hole it is just the extremal Reissner-Nordstrom (RN) black 
hole solution, as expected. A similar integral to (22) for the magnetic field reveals also that 
Pi stands for the magnetic charges. From the balancing gravitational and electromagnetic 
force the electric / magnetic charge are proportional to mass in accordance with 

Qi = rriie, (23) 

P . = m lV /r^2 ; (24) 

so that 

Qi + P? = < (25) 



III. GEODESIC EQUATION OF AN ELECTRICALLY CHARGED TEST PARTI- 
CLE 



In this section we seek for a solution to the geodesic equations of a test charge inside the 
field of a single static black hole located at the origin and for simplicity we shall assume 
A = 0, C = 1. The line element is given by 



1 



where 



n = i + 



m 



(26) 



(27) 



and m = \P = Q/e. The Lagrangian for a test particle with electric charge q and unit mass 
is 



t 2 n 2 



2fi 2 + 2 



•2,2/ q2 



+ ^t + qPcos6 0, 



(28) 
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in which a 'dot' stands for derivative with respect to the proper time r. This Lagrangian 
implies the following equations, and first integrals 

-i? + n = Q '°- (29) 

Q 2 r 2 sin 2 9 <p + qP cos 9 = (3 , 

where and (3q are two integration constants related to energy, angular momentum and 
the constant (5 is defined by — £p. 

We start with the 9 equation, by setting 9 = 9 . This leads to two different cases: 



A. (3 = cos# , (0 < < f) 

By taking /3 = cos#o and 9 = 9® one easily finds # equation is satisfied and <fi equation 
requires either 9 = or = 0. Here we exclude the case of 6q = and accept = 0. The r 
equation reduces to 

r + §r 2 = -^(ge-« fi). (30) 
The latter equation yields the following non-linear differential equations. 



1 . Case of ao = 

A specific analytical solution can be found by setting ao = 0. This choice leads to 

Q' 

r + —r 2 = (31) 

u L 

which reveals 

r (r) = m LambertW (e Ar+B ) (32) 

in terms of the LambertW(.) function jsj. The constants A and B can be fixed so that 
r (r) = is reached in a finite proper time. The ordinary time t is also expressed in terms 
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of the r by 

t(r) = % ( Lambert W (e Ar+B ) - - j— -. p , + 2 ln(LambertW (e Ar+B )) ] + C 

A\ v ; Lambert W (e Ar+B ) v " ) 

(33) 

in which C is another integration constant. In terms of the coordinate time t, x(t) satisfies 
the differential equation 

^ 1+I )$- 2 (I) 2 = ° < 34 > 

which can be studied numerically. 

2. Case of a / 

For the general case of a ^ 0, we introduce the new parameters 

r = mx, qe = m (B + A) , a — mA, (35) 

into the Eq. (30) to get 

x + %x 2 = -%A (B + A- Att) ->■ 

x(x + l) 2 x-(x + l)x 2 =A(Bx-A). (36) 

With particular boundary conditions we plot x (r) in Fig. la directly from this differential 
equation. Now, in order to obtain a particular solution we introduce the ansatz (from the 
analogy of a velocity dependent potential) 

oo 

x 2 = £a fc x fc (37) 

k=0 

which leads to 

oo Y 



x=J2^a k x k -\ (38) 
By substitution into (36) we get 

a = A 2 , <n = -2A(A + B), a k = - (2a k -i + a k - 2 ) = (-l) fe+1 [(k — 1) a + kai] , 

(39) 

which imply 



oo oo oo 

i 2 = E (-1) i(k - 1) a + k ai ] x k = a E (-1) + (* - 1) ^ + «i E ("1) kx k 

k=0 k=0 k=0 



(40) 



a I 1 - ( T— — ) I + °1 



x \ [(1 + 2x) ao + a;ai] 
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This easily gives 



and therefore 



±1 



x 



1 + x 



■y/A(A-2Bx) 



T + C 1 = 



±1 



3AB 2 
(Ci = constant). 



(A + B (3 + x)) v / A(A-2Bx) 



(41) 

(42) 
(43) 



This gives the relation between the proper time and the position of particle for any value of 



a . Next, by using the t component of the geodesic equation we find 

, / dr\ 2 d 2 r Q' 

{2qe - 3a ty ft' ( — J + Q (qe - a Q) — = -«o^, 

where Q is still given by (27). Here also we rescale our variables as 

r = mx, qe = m (B + A) , a — mA, t = ml 

to get 



72 

x (1 + xf (Bx -A)-^-(l + xf ((2B -A)x- 3A) 

LtL 



dx 
dt 



Ax 4 



which has the exact solution 

1 



±t + C 2 



3VA 



— 3 2A 
VA - 2Bx { - -6-x + — )+21n 
x B 



VA+VA- 2Bx 



VA-VA- 2Bx 



(44) 



(45) 



(46) 



(47) 



(A^O), (A>2Bx), (C 2 = constant).. 



We can easily observe that for x — > 0, t — > oo as expected for a distant observer; Fig. lb 
reveals this fact. 



3. The case of pure magnetic charge 

The case of pure magnetic charge can be obtained by setting e = 0, or equivalently 
B = —A. This leads to 

a = A 2 , ai = 0, (48) 

and therefore 

±\A\ 



X 



1 + x 



VT+2x, (49) 



t + C 3 = ^(2 + x)VT+2x~ (50) 
(C 3 = constant). (51) 
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The latter equation leads to 



x(t) = ~^/3^WTT^- T7 = L= = ) (52) 



in which 

a = ±\A\r + C 3 . (53) 



Eq. (45) becomes now 

— r M • rY' 

dt 2 V dt 



- X (l + x ) 5 ^ + 3(l + x) 4 (^] 2 = x 4 (54) 



with exact solution 



±t + C 4 = -Vl + 2x ( - -x-8) +21n 

3 \ x 



Vl + 2x + l 



(55) 



y/l + 2x-l 

(C 4 = constant). (56) 
^. T/ie case of pure electrically charged black hole 

By choosing e = 1, P = in the Lagrangian (28) we obtain a reduced set of geodesies 
equations. The 9 = case implies automatically that — | and = = /3q. This is nothing 
but same as (30) with the additional condition of e = 1, and the resulting geodesies motion 
obtained above. Thus, in class-A geodesies, P = case doesn't show a significant difference 
from the P^O case. 



B. = l/cos# , (0 < #o < f) 

After setting 9 = 9 , in order to solve 9 equation one can also choose (3 = 1/ cos9 . This 
choice in <p equation leads 

0= oS^' (57) 

and r = mx equation reads 

r + -r 2 = (qe - a Q) + ta n 2 9 . (58) 

This choice does not change the t equation. Now we use the same change of variables (35) 
together with 

q = mq,P = mP, (qP^J 2 tan 2 9 = C 2 (59) 
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under which the r equation takes the form 

where Q = 1 + -, and Q' = d x Q. It is observed that the last term on the right hand side is a 
direct contribution of the magnetic charge with marked distinction from the pure electrically 
charged black hole case. We note that by some manipulation on the equation, one gets 

• qP 1 V 



cos^ofi 2 ^ 2 Q?x 2 



(61) 



where T> = ^ . By transforming the independent variable from the proper time r to the 
azimuthal angle the orbit equation takes the form 



x 



or equivalently 



d<t> 2 x(l + x) V # / (l + x)£> 2V v /y D 2 (l + x)' 

Fig. 2 gives a numerical plot of x ((f)) , under the boundary conditions x (<p)\, =0 = 1 and 



4>=o 

dx((j>) _ q 

=0 

C. Generalization to two-centre black holes 

In this section we try to extend the result found for single black hole to double-black hole 
system. To do so we consider two identical black holes at (0, 0, h) and (0, 0, —h) , and the 
test particle is placed at a distance, far from the black holes such that one can write the 
metric function, up to the third order, as 

n=1 + 2m + mftW*-l) +0 /AV 
r r A \r J 

It is easily seen that for h — > the metric goes to the extremal RN black hole with both 
electric and magnetic charges and total mass 2m, as it should. The Lagrangian of the system, 
up to the same order of approximation, from the potential (2) can be written as 



i 2 n 2 



2il 2 



•2 i 2 



(^ 2 + sin 2 ^ 2 )]+^t+gP 



3h 2 sm 2 6\ „ [h 
oosfl ! 2 — ) +0 



0. 
(65) 
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This leads to the following geodesic equations (with integration constant ao and (3o) 

i= (qe- a tt) Q, (66) 



ft V sin 2 9 (f) + qP 



3h 2 sin 2 9\ „ (h 

cos6 ( 2 + o - 

r z I \r 



(67) 



f + 2 7T 



[qe - a iZ) ^ + 



r 2 + r 2 (fl 2 + sin 2 £0 2 ) 



rf^ + sinW) + ^ 



a 2 n 

.6h 2 sin 2 9\ //i 4 

COS 9 ( l+O -r- 



+ 



r 2 + r 2 ^ + sin 2 ^ 



0, 



+ ft V cos sin 6> 



(68) 



<?eft e 
ft 2 



i + gP 



sin0(2 + ^(3cos 2 0-l)') + o(* 



0. 



From (63-66) it follows that 



• (go-2gPcos0) Q (^h 
ft 2 r 2 sin 2 fl \ r 



— (r 2 Q 2 e) =Q 

(qP? / „ 3/i 2 sin^cos2# ... 

' /3 cos - 2 ^ /3 - 2 cos + 



Q2 r 2 s j n 3 £ 



qe 

qe - a Q — — (qe - a ^) + ^ 

,2 - 2 , 



f 2 + r 2 (fl 2 + sin 2 9 2 
3/i 2 cos sin 2 9 



(69) 



(70) 



+ 



(71) 



^3 1 ^ l / 



2m 3m/i 2 (3 cos 2 - 1) ^ //i 4 

= S- ^ + O — 



fi^-6 m/l2si ^ CQSg +0^) , (72) 

(73) 

so that the latter expressions satisfy the integrability condition, Q dr = ft r6 » within the range 
of approximation. We choose now, similar to the first case of single black hole case, the 
particular angles 



7T 



9 = ->■ 9 = -, 
2 

6 — ->■ 6 — 6o, 



(74) 



which give 



„ 2m 3mh 2 
O ~ i 



(75) 



ft ~ 1 + 



2m m/i 2 
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We see that the and 9 parts of the equations are trivially satisfied (by considering the 
approximation up to the third order) and the two remaining equations, i.e., r and t parts 
reduce to the same set of differential equations which were solved in the previous section, 
i.e., 

4 + l a a <" (76) 

fl o = 

r + -r )^-(qe-a n)—. 

It should be noted also that here the problem yields a different solution, because the metric 
function is different. Another special choice of interest to be considered here is given by 

e = o ->■ e = o, (77) 

= 0, 



which implies 



n e ~ 0, (78) 
^ 2m 6mh 2 

\l r ~ 

^ 2m 2mh 2 

+ — + — — . 

- + f = Q '°- (79) 

2qP = /?„, 



n J ^ 1 a 2 

These equations also make almost same set of equations as before. Let us add that gen- 
eralization to multi-coaxial black hole case (say, along the z-axis) can be treated more 
appropriately in the cylindrical polar coordinates. In these coordinates the electro magnetic 
potential ansatz takes the form 



with 



A = ^dt + Y. 1>A - l ' ] =<lo (80) 
1 ^ P 2 + (z- Zl f 



^ = i + E^= m (81) 

\2 
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and the constraint reads as before, namely 

2 1 

e + p = 

777,3 = 

This describes an infinite array of MP black 
magnetic charges, and the line element is gh 



|AP 4 | = 5i. 

e 

holes, each at z — Zi, with coupled electric and 
^en by (1). 



IV. CONCLUSION 



We extend the electrically charged MP black holes to the dyonic case which possesses 
both electric (Qi) and magnetic (Pi) charges. Superposition principle provides us multi- 
black holes where the mass (mi) of each black hole satisfies Pf + Q\ — mf. The charges 
are scaled by a parameter e (0 < e < 1) which regulates the effective charges of both types. 
Under such restriction only we were able to obtain such dyonic solutions. In order to find the 
interior charge content of the black hole we provide a detailed analysis of geodesies. Exact 
particular integrals are available in some cases but for the general treatment we resort to 
the numerical integration and two-dimensional plots. The orbit equation reveals also the 
hovering of a test charge around a dyonic black hole. By a detailed analysis it seems possible 
that we may identify the charge constituent of a MP black hole. In a more heuristic argument 
a magnetically charged black hole may be identified as a magnetic monopole which, so far 
has not been detected in our observable universe. As a final remark we wish to add that 
with the inclusion of time in the metric a' la [4| collision problem of magnetic MP black 
holes can be investigated. 
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Appendix: 

Regarding the MP line element, the non-zero energy momentum tensor and Einstein's 
tensor components are 



Tl = - {Fl + F 2 y + F 2 Z ) - — [(F% + F 2 XZ + F 2 Z )] , (la) 
T x - + F 2 + F 2 ) + — (F 2 +F 2 - F 2 ) 

1 x — \ r tx ' r ty ' r tz) ' ^4 \ r xy ' r xz r yz) i 

T y = (p2 _ p2 +p 2\ , J_ ip2 _ p 2 +F 2\ 
± y \ r tx r ty * r tz) * ^4 \ r xy r xz ' r yz) i 

T z = (F tx + F ty — F tz ) + — (—F xy + F xz + F yz ^j , 

2 

Tx — T y = —2F tx F ty + — (F xz F yz ) , 

2 

T x — T z = —2F tx F tz + (F xy F zy ) , 
T; = Tj> = -2F ty F tz + A (F yx F zx ) . 
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^ = -[f + ^) » (lc) 

rpx _ (2 , i\ (-ni + & y + ni) 
x V a 2 / ^ 4 



A 2 / ft 4 



Tl = ( + - 
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A 2 / ft 4 



= t; = - ( 6 2 + 
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Figure captions: 

Fig la: Freely falling charged particle into the dyonic black hole as a function of proper 
time. In a finite proper time the particle reaches the horizon, as expected. With the magnetic 
charge on the black hole, the test particle plunges into the black hole in a shorter proper 
time. The infall gets delayed for a weaker magnetic charge. 

Fig lb: The free fall motion of a test charge is observed from a far distance. It takes an 
infinite coordinate time to reach the horizon and the magnetic charge has little effect in the 
process. 

Fig 2: The oscillatory motion of a test charge around a dyonic black hole, x ((/>)(— r ((f))) is 
plotted versus the azimuthal angle. Our boundary conditions are such that x (0 = 0) = 1 and 

4f =0, the rest is determined by the differential equation of orbit. An exact solution, 

d<p 0=0 

which is not at our disposal, should definitely reveal much more than our numerical analysis. 
The effect of the magnetic charge on the behavior of the test particle is evidently visible. 
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This figure "Figla.jpg" is available in "jpg" format from: 
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This figure "Figlb.jpg" is available in "jpg" format from: 
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